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Abstract 

Making use of a Rice-like series expansion, for a class of stationary Gaussian processes the asymp- 
totic behavior of the first passage time probability density function through certain time-varying 
boundaries, including periodic boundaries, is determined. Sufficient conditions are then given such 
that the density asymptotically exhibits an exponential behavior when the boundary is either asymp- 
totically constant or asymptotically periodic. 
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1 Introduction 

First-passage-time (FPT) probability density functions (pdf 's) through generally time-dependent bound- 
aries play an essential role in many applied fields including the stochastic description of the behavior of 
various biological systems (see, for instance, [8], [21], [22], [25], [28] and the references therein). Investi- 
gations have essentially proceeded along the following three main directions: (i) to search for closed-form 
solutions under suitable assumptions on the considered stochastic processes and on the boundaries (see, 
for instance, [7], [11], [16], [18], [20]); (ii) to devise numerical algorithms to evaluate FPT densities (see, 
for instance, [2], [3], [4], [5], [6], [13], [14], [15], [17], [29]) and (Hi) to analyze the asymptotic behavior of 
the FPT densities as boundaries or time grow larger (see, for instance, [19], [23], [24], [30], [31]). The 
present paper, that falls within category (Hi), is the natural extension of previous investigations carried 
out by us for the class of one-dimensional diffusion processes admitting steady state densities in the pres- 
ence of single asymptotically constant boundaries or of single asymptotically periodic boundaries ([19], 
[23], [24]). In such cases, computational as well as analytical results have indicated that the FPT pdf 
through an asymptotically periodic boundary is susceptible of an excellent non-homogeneous exponential 
approximation for large times and for large boundaries. 

However, if one deals with problems involving processes characterized by memory effects, or evolving 
on a time scale which is comparable with that of measurements or observations, the customarily assumed 
strong Markov property does not hold any longer; hence facing FPT problems for non Markovian pro- 
cesses becomes unavoidable. As is well known, for such processes no manageable equation holds for the 
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conditional FPT pdf: only an excessively cumbersome series expansion is available when the process is 
Gaussian, stationary and mean square differentiable (cf. [26], [27] and the references therein). 

Due to the outrageous complexity exhibited by the numerical evaluation of the involved partial sums 
on accounts of the analytical form of the involved terms, a totally different approach has been recently 
undertaken in order to obtain information on the asymptotic behavior of the FPT densities for a class 
of normal processes. This consists of a simulation procedure (see [9] and [10]) implemented to generate 
sample paths and to estimate the corresponding FPT densities. Extensive computations have thus been 
performed to gain some insight on the behavior of the FPT pdf through varying boundaries. The results 
of the simulations, obtained by means of a parallel supercomputer CRAY T3E, have indicated that for 
certain periodic boundaries not very distant from the initial value of the process, the simulated FPT pdf, 
g(t), soon exhibits damped oscillations having the same period of the boundary. Indeed, to a high degree 
of accuracy, g(t) can be represented in the form 

~g{t)~ht)e-«\ (1.1) 

with a and (3(t) specified by means of the data obtained via the performed simulations [12]. Note that 
(1.1) can be thrown in the equivalent form 

g(t) ~ a(t) ex p{~^ a(r)dr}, (1.2) 

where a(t) > is a periodic function having the same period of the boundary. Hence, for periodic 
boundaries, even though not very distant from the initial position of the process, the estimated FPT pdf 
appears to admit a non-homogeneous exponential approximation. 

In the present paper the relevance and the validity of such an unexpected numerical result is confirmed. 
Indeed, it will be proved analytically that the non-homogeneous exponential approximation (1.2) holds 
for a wide class of stationary Gaussian processes in the presence of boundaries that either possess a 
horizontal asymptote or are asymptotically periodic. 

In Section 2 we shall briefly recall some basic notation that will be used throughout this paper; in 
Section 3 we shall assume that the boundaries possess a horizontal asymptote, and in Section 4 that 
they are asymptotically periodic. Finally, in Section 5 for a stationary Gaussian process with zero mean 
and damped oscillatory covariance, the simulated FPT pdf g(t) is compared with the non-homogeneous 
exponential approximation for the FPT pdf. 



2 Mathematical background 

Let {X(t), t > 0} be a one-dimensional, non-singular stationary Gaussian process with mean E[X(t)] = 
and covariance E[X(t)X(r)} = j(t - r) such that 7(0) = 1,7(0) = and 7(0) < 0. Then X(t), the 
derivative of X{t) with respect to t, exists in the mean-square sense. Let S(t) € C 1 [0, +00) be an 
arbitrary function such that -X'(O) = x < S(0). Then, 

T=M{t:X{t)>S(t)}, X{0)=x o (2.1) 
is the FPT random variable and 

g(t\x ) = ^P(T < t) (2.2) 

is the FPT pdf of X(t) through S(t) conditional upon X(0) = x . For all n e N and = t < h < . . . < t n 
we denote by W n (t\, . . . , t n \x ) dt\- ■ ■ dt n the probability that X(t) crosses S(t) from below in the intervals 
(t\,ti + dti), . . . , (t n ,t n + dt n ) given that X(0) = x - As shown in [26], the functions W n can be expressed 
as 



r + ao 


r+00 


r+00 


1 d vi 


/ dy 2 -- 




/S(ti) 


fs(t 2 ) 


Js(t n ) 



n i# - 



i=l 



xp2n[S(ti),t 1 ; . . . ;S(t n ),t n ;yi,ti, . . .;y n ,t n \x ] dy n , (2.3) 
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where P2n{xi,t\; ■ ■ ■ , x n , t n ; . . . ; y n , t n \xo) is the joint pdf of 2n random variables X(t\), . . . , X(t n ), 

Y(h) = X(h), . . ., Y(t n ) = X(t n ) conditional upon X(0) = x : 



P2n(xi,ti; . . .;x n ,t n ;yi,ti; . . .;y n ,t n \x ) 

2n 



X cxp 



2 |A27i+i(ti, . . . , t n )\ 



(27r)™|A 2 „ +1 (t 1 ,...,i„)| 1/2 
^2 k+i,j+i(h, ■ ■ ■ ,t n ) [xi - x j(U)] [xj-x i(tjj\ [. 

i,.7 = l ' 



(2.4) 



Here h+i,j+i(ti, . . . , t n ) denotes the cofactor of the element \ i+1 j + i(ti, . . . ,t n ) of the covariance matrix 

A2n+l(tl, ■ ■ ■ , t n ) Of X(0), X(h), X (t n ) , X (h) , ...,X (t„), i.e. 

f £[A(^)A(t,)] = 7(tj - tj) = y(tj - U), 

(i = 0,1,..., n,j = 0,1,..., n) 



K+i.j+i(ti, ■ ■ ■ , t n ) — < 



E[X(ti)X(tj- n )] = -7(ti - *,•_„) - 7fe-n - 

(i = 0, 1,.. .,n, j = n + 1,. ..,2n) 

E[X(t^ n )X{t 3 )} = 7 (t 4 _„ - tj) = -7(tj - ti_ n ), 

(j = n + 1, . . . ,2n,j = 0, . . . ,n) 

E[X(ti- n )X(tj-„)] = -j(ti- n - tj- n ) = —j(tj- n - U- n ), 

(i = n + 1 . . . , 2n, j = n + 1 . . . , 2n) , 



(2.5) 



Xj, .- 



7(*i) 



Di- n , i = n + l,n + 2,...,2n 



-{ 



l(U), i = l,2,...,n 
jiU-n), i = n+l,n + 2,...,2n 



and \A\ denotes the determinant of a matrix A. Substituting (2.4) in (2.3) one has: 



W n (ti, . . .,t n \x Q ) = 



(27r) n |A 2n+ i(ii,...,t r . 



,1/2 



r+oo i> 



d^ 2 



(2.6) 



/*+oo n 

■ [T [& ~ ^(**)] • • • .*n;£l, • • • ,€n)K(tl, ...,t n ;£i,.. .,£„) d£ n , 

U(t n ) i= i 



where 



^(*i) = S(*i)-a:o7(ti) (* = l,2,...,n) 



. . . ,* n ;£l, • • • ,£n) = cx p{~Tu 77 77TT '*+n+l,j+n+l(<l) • • • , *n)& £j } 

2 |A 2n +l t„j| Z*-^ J 



^(ii, . . . ,t n ;£i) • • ■ , £, 



n) = expj-- 



(2.7) 



2 |A 2n+ i(ti, . . . , t n )\ I 

n n -in 

+ ^ l i+l,3+n+l(tl,---,t n )'4>(ti)£j+ ^ l i+n+l,j+l{tl,---,t n )lp(tj)£i >. 
i,j = l »,j=l 

As shown in [26], g(t\xo) can be expressed as the following Rice-like series: 

oo „t ft ft 

g(t\x ) = W 1 (t\x ) + J^i-iy / dti / dt 2 --- / . . . ,*i,tjas ) dti , (2.8) 

i=1 JO Jti 
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with x < S(0). 

Let now a r denote the partial sum of order r of the series in (2.8). Then, for each t > the partial 
sums of even order give a lower bound to g, whereas the partial sums of odd order provide an upper 
bound to g. Since the evaluation of the partial sums is very cumbersome because of the complexity of the 
functions W n and of their integrals, a first approximation of FPT density can be carried out by evaluating 
Wi(t\x ). The explicit expression of Wi(t|a;o) (cf. [27]) is: 

Wi{tM -2n[i-^( t )] exp v w^m\ / ( } 

f a 2 (t\x )\ pK { <r(t\x ) 

cx p i s — r - \/ o Erfc 



where 



|As(t)| = -7(0) [1-7 2 W]-[7W] 2 

^M-f^#V /2 |^) + 7W[7(t) ^'" ] 



(2.10) 



|Aa(t)| J W 1-7 2 W 



and 



2 f + r ?r ~i 

Erfc(z) = -^= J exp{-^jdy, z e R. (2.11) 

We stress that although (2.8) gives a formal analytical expression for the FPT densities through 
arbitrary time-dependent boundaries, no reliable numerical evaluations appear to be feasible due to the 
complexity of (2.6) and (2.7). Furthermore, for all t > the first-order approximation Wi(t|a;o), that 
provides an upper bound to the FPT pdf in (2.8), is a good approximation of g only for small values of t. 

3 Asymptotically constant boundary 

In this Section we consider the FPT problem for an asymptotically constant boundary 

S(t) = S + Q(t), t>0, (3.1) 
with So € R and where g(t) € C 1 [0, +oo) is a bounded function independent of So and such that 

lim g(t) = and lim g(t) = 0. (3.2) 

l — > + oo t — >+oo 

The following proposition shows that under suitable hypotheses on the covariance j(t), the function 
Wi(t|a;o), given in (2.9), approaches a constant value as t increases. 



Proposition 3.1 If 



then 



lim 7(i) = and lim j(t) = 0, (3.3) 



R(S ) := lim W^xq) = ^ ^ exp (-^ } . (3.4) 



t^+oo 2 7T 2 
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Proof. Due to (3.2), from (3.1) it follows 



lim S(t) = S , 

t— >+oo 



lim S(t) = 0. 

i— >+oo 



Furthermore, by virtue of (3.3) and (3.5), from (2.10) one has 



lim |A 3 (i)| = -7(0), 



lim a(t\xo) = 0. 



(3.5) 



(3.6) 



Taking the limit as t — > +oo in (2.9), and making use of (3.5) and (3.6), we are finally led to (3.4). This 
completes the proof. ■ 



We note that R{S ) > for all S Q e R and 



lim R(S ) = 0. 



(3.7) 



The nonzero asymptotic value given by (3.4), per se indicates the inadequacy of Wi(t|a;o) to provide a 
valid approximation to g(t\xo) for large time. On the contrary, the goodness of such an approximation 
for small times is confirmed. 



Theorem 3.1 Let S(t), i > 0, be bounded and such that (3.1) and (3.2) hold. If 



lim 7 W = 0, lim 7(i)=0, lim 7 (t) = 



and 



lim 



R(S ) 



with R(Sq) defined in (3.4), then 



lim 



1 



So 



Xq 



x < 5(0). 



So^+oo R(S ) a \R(S ) 
Proof. Since R(Sq) > for all So € R, changing t in t/R(So) in (2.8) we obtain: 



i2(So) V^(^o) 



1 



£o = 



R(S ) \R(So) 



t 



x 



+ 



and hence, 



-. oc — £ 



dti 



R(S a ) 



di 2 ■ ■ ■ 



tl 



H(S ) 



R(So) 



xo dU 



1 



-.9 



R(So)*\R(S ) 



t 



Xo = 



R(So) \R(So) 



after the change of variables Tj = R(So) for i = 1,2,.... 
We now prove that for all n = 1, 2, . . . one has: 



)'•••' R(S )'R(S ) 



x dn 



lim 



So^ + oo [R(S )Y 



n 



XQ 



R(So) ' ' R(So) 

(0 < n < r 2 < . . . < r„). 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 
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To simplify the notation, we set: 



R(S ) , --- , " n R(S ) 
Hence, recalling (2.6), for n = 1, 2, . . . one has: 



(0 < 1?! < 2 < . . . < n ). 



(27r)"|A 2 „ +1 (0 1 ,...,0„)| 1/2 



(3.13) 



(3.14) 



/ #1 / ^•••/ nki-V'(^) ff(<?l,-"^n^l,...,Cn) 

./i/>(#i) Ji/>(#2) Jiif^l „•_-, L 

1 



[fl(So)]' 



i= i 

^ ■ ■ ■ j £l> ■ ■ ■ i £n) l dS,n 



where the functions ip, H and K are defined in (2.7). 

We shall now make use of the following relations (see Appendix I): 

lim V>(^)=0 (i = 1,2,..., n), 



lim 



So- 



A 2 „ +1 (tf 1 ,...,tf„) =[-7(0)]", 



lim 

So— >+oo 





(3.15) 




(3.16) 


7(0)] g^}' 


(3.17) 


(2tt)» 

[- 7 (0)]»/2 


(3.18) 



Taking the limit as So — > +oo in (3.14) and recalling (3.15), (3.16), (3.17) and (3.18) one then obtains: 

: W„(i?i,...,1? n | x ) (3.19) 



So^oo [iJ(S'o) 



{^2pK0)l} c 



[-7(0)] 



/•+oc 

J C exp 



Since 



/■+00 /• (-2 ~\ r + oo 



2 [-7(0)]. 

relation (3.12) immediately follows from (3.19). Due to (3.12), taking the limit as So — > +oo in (3.11), 
one finally obtains: 



lim 



1 



(3.20) 



5o^+oc i?(5 ) a \R{S ) 

that identifies with (3.10). The proof is thus complete. ■ 
The following corollary is an immediate consequence of Theorem 3.1. 
Corollary 3.1 Under the assumptions of Theorem 3.1, for So — > +oo one has: 

g(t\x ) ~ R(S ) exp{-i?(S ) t}, Vt > 0, (3.21) 

with R(S ) defined in (3.4). 

This Corollary expresses the asymptotic exponential trend of the FPT density as the boundary moves 
away from the process' starting point. 
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4 Asymptotically periodic boundary 

The FPT problem in the case of an asymptotically periodic boundary will be the object of the present 
Section. More specifically, we shall focus our attention on boundaries of the form 

S(t) = S + g(t), t>0, (4.1) 

where S eR and g(t) € C 1 [0, +oo) is a bounded function independent of So and such that 

lim g(t + kQ) = Z(t), lim g(t + kQ) = Z(t), (4.2) 

k— >oo k^oo 



where Z(t) is a periodic function of period Q > satisfying 

/o 



/ Z(t) oIt = 
Jo 



0. 



(4.3) 



Proposition 4.1 If (3.3) holds, then 



R[Z(t)} := lim Wl (t + kQ\xo)= VEM ex J J S * + Z ^ 

fc^oo Z IT \ Z 



exp 



2[-7(0)] / 

Proof. From (2.9) for k = 0, 1, ... we have: 

M/^tm ^ \A 3 (t + kQ)\^ 
Wl{t + kQlX0)= 2^WTkQ)] e * P 

f <7 2 (t + fcQ|x ) 

exp <^ 



Z(t) Erfc 



Z(t) 



2 [-7(0)] ' [-7(0)], 

f [S(t + fcQ) -.T 07 (t + fcQ)] 2 ] 

1 2[l- 7 2 (t + fcQ)] J 

. , n , x w f (<?{t + kQ | ) 
a(t + kQ\ x ) Erfc ' 



2 J V 2 " v ~ ' " """ """ V V2 

where a(t\xo) and |A 3 (t)| are defined in (2.10). By virtue of (4.2), from (4.1) it follows: 



lim S{t + kQ) = S + Z(t), 

k — *oo 



lim S{t + kQ) = Z(t). 

k— >oo 



Furthermore, due to (3.3) and (4.2), from (2.10) one has: 

lim |A 3 (i + kQ\x )\ = -7(0), lim a(t + kQ\x ) 



k— *oo 



z(t) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



Taking the limit as k — > +oo in (4.5), and making use of (4.6) and (4.7), we are then led to (4.4), which 
completes the proof. ■ 



Remark 4.1 For all t > the function R[Z(t)] defined in (4.4) is a positive, periodic function with 
period Q. Furthermore, R[Z(t)\ can also be written as 



R[Z{t)] = ^ = exp(- [ ^ + Z(t)]2 l r^-Zfi)] exp(- r ^ „ \ d£. 

1 UJ 27rypK0)] P \ 2 /y i(t) [? WJ P \ 2[-7(0)]/ 



(4.8) 
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Proof. Since Z(t) is a periodic function of period Q, for all k = 0, 1, . . . one has Z(t + kQ) = Z{t) and 
Z(t + kQ) = Z(t). Hence, from (4.4) it follows that R[Z(t + kQ)] = R[Z(t)] (k = 0, 1, . . .), i.e. R[Z{t)\ 
is a periodic function with period Q. 

We shall now prove that R[Z(t)) > for all t > 0. To show it, let us re-write R[Z(t)} as 



fl[Z(t) ] = ^Zpexp 



2tt 





r £(*) l 




l-V2[-7(0)]J 



with 



e~ x * dx. 

We note that A(y) > for all y < 0. Furthermore, if y > the following inequality (cf. [1]) holds: 



(4.9) 



(4.10) 



y + yjy 1 + 2 



< 



e x dx < 



Hence, one has 



\ y 2 + --y 



_ H 2 yj y 2 + 2 - y 

y + x/v 2 + 2 



(y > o), 



which again implies A(y) > 0. Therefore, for all y € R there holds A(y) > 0. Recalling (4.9), it 
immediately follows R[Z(t)] > for all t > 0. 
Finally, to prove (4.8) we note that 



r+oo 

Jz(t) 



£ - z(t) 



exp 



i 2 



2 [-7(0)] 



= [-7(0)] A 



Z(i) 



(4.11) 



with defined in (4.10). Hence, making use (4.11) in (4.9), equation (4.8) immediately follows, as it 
had to be proved. ■ 



Since Z(t) does not depend on So, from (4.4) it follows: 

lim R\Z(t)] = 0. 



(4.12) 



Proposition 4.2 Let 



1 f Q 

a = a(S ) := - / R[Z(r)} dr, 



(4.13) 



with R[Z(t)] defined in (4.4). Then, there exists a non-negative monotonically increasing function ip(t) 
which is a solution of 



/ R[z( 

Jo 



r)l dr = at, 



Vi > 



such that 



¥>(0) = 0, 
lim iy2(i) = +oo 

V>(t + fcQ) = v»(t) + (A; = 0,1,...). 



(4.14) 



(4.15) 
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Proof. In Remark 4.1 we have proved that R[Z(t)\ > 0, Vt > 0; hence, from (4.13) it follows a > 0, 
and from (4.14) one has <p(0) = and <p(t) > 0, Vt > 0. Let h(t) be any primitive function of R[Z(t)]. 
From (4.14) we have h{<p{t)\ = /i(0) +at. Since R[Z(t)\ > 0, Vt > 0, /i(i) possesses an inverse, and hence 
ip(t) — /i -1 [/i(0) + at}. Furthermore, since a > 0, from (4.14) one has 

i^'wbm^ (t>0) - (4 ' 16 » 

Therefore, ip(t) is a monotonically increasing function for all t > 0. Furthermore, since R[Z(t)\ is a 
positive function, the second of (4.15) holds. We now remark that from (4.14) one has: 



r<p(t+kQ) p<p(t) r-ip(t + kQ) 

/ R[Z(t)]oIt= R[Z(t)]oIt+ R[Z( T )}dT = a{t + kQ), (4.17) 

JO JO Jtp(t) 



rip{t+kQ) r <p(t) r<p(t+kQ) 

Mi) 

or, due to (4.13), 

r <p(t+kQ) r Q r kQ 

M*) 



p<p(t+kQ) r Q rkQ 

/ R[Z(T)]dT = kaQ = k R[Z(t)}oIt= R[Z(t)} dr, 

J<p(t) Jo Jo 



(4.18) 



where the last equality follows since R[Z(t)\ is a periodic function with period Q. Relation (4.18) finally 
implies the last of (4.15). The proof is now complete. ■ 



Proposition 4.3 For all t > one has 

(o K£) >0 ' 



R 



(Hi) lim (p 



+oo, 



(iv) lim \ip(— ip(—\\ = +oo (0<T<t). 



Proof. Since <^(t) is a non-negative function and a > 0, condition (%) follows from (4.13), while from 
(4.16) immediately one obtains (ii). Making use of (4.12), from (4.13) we have 



lim a = i lim [ Q R[Z(r)} dr = 0, 



(4.19) 



that, due to the second of (4.15), implies (Hi). Finally, making use of the mean theorem of Calculus one 
has: 



\aJ \aJ a \ 



t - T 



R 



(o < r < i < t), 



(4.20) 



where the last identity follows from (ii). We note that, due to (4.12), there holds: 

t 



lim R 

So— >+oo 



o. 



(4.21) 



so that (iv) follows after taking the limit as So — » +oo in (4.20). The proof is thus complete. 
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The following theorem then holds. 

Theorem 4.1 Let S(t),t > be given in (4.1) with g(t) G C^O, +oo) a bounded function such that 
(4.2) and (4.3) hold. If the covariancc function -f(t) satisfies (3.8) and if 



lim 



S + Z[ v [ - 



o, 



(4.22) 



with ip(t) defined in (4.14), then 



lim 

So— >+oo 













fV-Y 

dt \a J 
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K=) 




= e 









x < 5(0). 



(4.23) 



Proof. From (i) and (in) of Proposition 4.3 it follows that ip(t/a) can be viewed as a scaled time. 
Changing t to <p(t/a) in (2.8), we then obtain: 



d ( t 









'd 


a/ 






dt 


U) 






-) 



dti 



■ V 



dt 2 ■ ■ ■ 



x 



dt 



W i+1 [ t 1 ,---,t i ,ip[ — 



x dU . 



Hence: 



dt 



+ 



x 



dt 



W 1 



i=l J ° 



d ( T\ 



JTi-1 



dn \ a 



03 — 

dr\ V a 



dr\ 



f 



x 



d ( T2 



W i+ i a? — >•••,¥> - ],<P - 



Q: 



■03 — 

dr 2 \ a 



a 



dr 2 - ■ ■ 
x I dn , 



(4.24) 



after having performed the change of variables Ti — ip y — ). Due to (ii) of Proposition 4.3, (4.24) can 
also be written as 
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K=) 


x 


i? 




))] 



+ / dr i / dr 2-" ( 4 - 25 ) 

j=1 JO 7 T1 



J Ti-l 



w i+1 v 



xq 



R 



■R 



R 



dn . 



Let us now prove that for n = 1,2,... there holds: 



lim 

So^+oo 



W n [<p[ — 



1 a 



.T() 



1? 



V a 



I? 



Z\<p 



1 (0 < n < r 2 < . . . < t„). 



(4.26) 
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For simplicity of notation, we set: 



^ = K a) ' ' ' ' = ^(5) (0 < ^1 < ^2 < • • • < t?„). 
From (4.8), for n = 1, 2, . . . it then follows 



R 



» + oo ;>+oo /- + oo 

x / d£i / <*&■■■/ 

'z(i?i) iz(« 2 ) Jz(#n)7= 



(2 7r)"[- 7 (0)] n/2 



Hence, 



(ft)] ■■■i2[Z(0„) 



(2^)"[-7(0)]" /2 



exp< 


[i£[5o + ^)] 2 ' 


> W„ (<?i,- ■ 


■,ft 




/ #1 / d£ 

JZ(#i) iz(« 2 ) 


2-/. n 

»=1 


- ^ (<?<)" 


exp< 


I 2[- 


1 ™ 
-7(0)] S^J 





where, due to (2.6), one has 

W„(ft,---,1?„ I so) = 

/■+00 n 



., , d£i / d^ 2 ■ ■ ■ 

(2 tt)™ |a 2 „ +1 (ft, ... , i?„)| i/2 y^(*o y*(tf 2 ) 



with V', if and if defined in (2.7). 
We now note that (see Appendix II) 



lim 

So — >+oo 



A 2n+1 (ft,...,ft) =H(o)] n , 

lim^ H (ft, . . . ,0 n ;£i, . . . ,£„) = cxp|-^-L_ g£ j 
cxp|^[^o + ^(^)] 2 }a-(^,...,,?„; Ci,. ..,£„) 



lim 

So^+oc 



1. 



(4.27) 



|-^X>o + ^)] 2 } (4.28) 



(4.29) 



(4.30) 



Taking the limit as S -> +00 in (4.29) and recalling (4.31), (4.32) and (4.33), for all n = 1, 2, . . 

W n (ft,- ••,!?„ I So) 
S ™ooi?[Z (ft)] (T? n )] 



(4.31) 
(4.32) 

(4.33) 
one has 
(4.34) 



lim 

So^+00 



r + OO 

/ #1 


r + oo 


r+oo n 

1. n 

V>(i?„) »=1 




exp< 


I 2[-7(0)]£^J 




r+oo 

^ 1 


f 00 d£ 2 ■ ■ ■ 1 

Z(# 2 ) J 


r+°° n 

n 

Z(i? n ) i=l 




exp< 


I 2 [-7(0)] 





t -TT U I ^0) 

=^11 vm 

i=i 



(0<ft <ft < ••• <ft), 
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where we have set: 



U(0i\x o ) 



exp 



2 [-7(0)] 



(4.35) 



f (*) 



+ 00 



exp 



2 [-7(0)] J 



We now note that from (4.10), (4.11) and the first of (2.7) it follows 

[S(0i)-a;o7(0i)] 2 



C (0< | ar ) = -7(0) exp 



2 [-7(0)] 



tt [-7(0)] 



5(0i) -107^) Erfc 



(4.36) 



^W = -7(0)exp{-^||}- 

Furthermore, recalling (3.8) and (4.6), one has 

U(di\x ) 



7T [-7(0)] 



Z(i?i) Erfc 



£(0i) 



V2[- 7 (0)] 



lim 



= 1 (* = l,2,...,n) 



so that (4.26) immediately follows from (4.34). Due to (4.26), taking the limit as So — > +oo in (4.24), 
one obtains: 



lim 

S ->+°o 
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K=) 











00 /*i /-i /-t 

= l + 2(-l)* / dn / dr 2 -- - / dr,. 



(4.37) 



Equation (4.23) finally follows from (4.37). The proof is thus complete. 
Corollary 4.1 Under the assumption of Theorem 4.1, for So — > +oo one has: 



p (t | a: ) ~ R[Z(t)} expj-^ fl[Z(r)] drj, Vt > 0, 



(4.38) 



with R[Z(t)] defined in (4.4). 

Proof. From fw) of Proposition 4.3 and from (4.23) we obtain: 



lim 



R 



x 



x < 5(0). 



(4.39) 



Hence, using (4.14), we have asymptotically: 



r 1 - 



X 



R 



Z[ V [- 



R 



r f<p(t/a) | 
exp|-y i?[Z(T)]rfr|. 



The asymptotic formula (4.38) then follows. The asymptotic non-homogeneous exponential behavior of 
the FPT density has thus been proved. ■ 
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Note that (4.38) can also be written as 

g(t\x )~ /3(t)e~ at , 
where (3(t) is a periodic function of period Q given by 



(3(t) =R[Z(t)} expjat-^ R[Z(t)} drj, 



(4.40) 



(4.41) 



with a defined in (4.13). Indeed, since R[Z(t)] is a periodic function of period Q, due to (4.13) and (4.14), 
one has: 

{r-t+nQ 
a (t + nQ) -J R[Z(T)}dTj (4.42) 

r rt+nQ ~| 

= 0{t) expja nQ - J R[Z(t)] drj 

J R[Z(T)]dT-J R[Z(t)} drj =P(t). 
Notice that in (4.40) is an asymptotic expression of g(t | xq) of the same form as (1.1). 
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Figure 1: Plots of the function Wi(t|0) and of the simulated FPT density g(t) for a zero-mean stationary 
Gaussian process originating at x n = having covariance (5.1) with a = w = 1 and for the constant 
boundary S(t) = 2. 



5 Analysis of a special case 

The purpose of this Section is to analyze the behavior of the FPT pdf for a stationary Gaussian process 
of concrete interest for certain applications. 
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Figure 2: For the boundaries S(t) = So = 2 and S(t) = Sq = 2.5, the asymptotic exponential approx- 
imation R(Sq) exp{— R(Sq) t] and the simulated FPT density g(t) are plotted for the Gaussian process 
of Figure 1. 

Let {X(t),t> 0} be the stationary Gaussian process originating at xo = 0, with zero mean and 
damped oscillatory covariance [32]: 



7 (t) = e 



- p -a\t\ 



COS I 



(ut) + — sin(w • \t\) 



(t e R), 



(5.1) 



where a and to are positive real numbers. Functions of form (5.1) can often be conveniently used to 
approximate experimental covariance functions that, starting from a unit initial maximum amplitude, 
asymptotically tend to zero with an exponential envelope. From (5.1) we see that 7(0) = 1. Furthermore, 
7(0) = and 7(0) = - (a 2 + lo 2 ) < since for t > there holds: 



7(t) - -i(-t) 



-a t 



(a 2 + lu 2 ) sin(a;i), 



7(t) = j(—t) = (a 2 + uj 2 ) a sin(u! t) — lu cos(cj t) 

to 

We finally note that lim j(t) = 0, lim j(t) = and lim j(t) = 0. 



(5.2) 



5.1 Constant boundary 

For the constant boundary S(t) = So, from (3.4) one has 



R(S ) = 



2tt 



,S 2 o/2 



(5.3) 



All forthcoming figures refer to the case a = to = 1 in (5.1). Figure 1 shows the plots of VFi(t|xo) given 
by (2.9) and of the simulated FPT density g(t) as functions of t for the constant boundary S(t) = 2. 
The first-order approximation Wi(t|a;o) is seen to provide an upper bound to the FPT pdf in (2.8), 
being a good approximation of g only for small values of t. From (3.4) it follows that, as t increases, 
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Figure 3: Same as in Figure 1 for the periodic boundary S(t) = 2 + 0.5 sin(27ri/3). 



Wi(t\x ) approaches the constant value R(2) = e~ 2 /(7T\/2) = 0.0304611. Making use of (3.21), we 
expect that R(So) exp{— R(So) t} provides a good approximation of the simulated FPT density g(t) as 
5o — > +oo. Indeed, Figure 2 shows the plots of the simulated FPT density g(t) for the constant boundary 
S(t) = Sq = 2 and of the function R(2) exp{— R(2) t}. Figure 2 also shows the plots of the simulated 
FPT density g{t) for the constant boundary S(t) = Sq = 2.5 and the function i?(2.5) exp{— R(2.5) t} with 



R(2.5) 



-3.125 



/(7r\/2) = 0.00988928. We note that starting from rather small times, g(t) is susceptible 



of an excellent exponential approximation already for positive boundaries of the order of a couple of units. 



5.2 Periodic boundary 

For the same stationary Gaussian process, we now assume that S(t) = So + B sm(2irt/Q). In this case 
Z(t) = g(t) = B sin(2 7rt/Q) and condition (4.3) is satisfied. From (4.4) one has: 



R[Z{t)] = 



2n 



■cxp< - 



s 2 (t) 



cxp< - 



2(a 2 +uj 2 ) 



2(a 2 +Lu 2 ) 



Z(t) Erfc 



Z(t) 



v/2 (a 2 +iu 2 ) 



(5.4) 



From Remark 4.1 it follows that the function (5.4) is positive and periodic with period Q. 

Figure 3 shows the plots of VFi(£|0) and of the simulated FPT density g(t) for the periodic boundary 
S(t) = 2 + 0.5 sin(27r£/3). We note again that the first-order approximation W\ provides an upper 
bound to the FPT pdf given by (2.8), though being a good approximation of g only for small values of t. 
Furthermore, as k increases W\{t + kQ \ 0) tends to the function (5.4), i.e. as t increases Wi(f|0) becomes 
periodic with the same period of the boundary. Since 



lim 



So + Z[if[ - 



= lim 



B sin^2 7r 



<P(t/a) 



S + Bsin(2n^M) 



0, 



Q 



condition (4.22) is satisfied, and thus the asymptotic formula (4.38) holds. 
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5 10 15 20 25 30 35 40 

Figure 4: For the same process as in Figure 1 and for S(t) = 2 + 0.1 sin(27r t/3), the function 
R[Z(t)\ exp{— J Q R[Z(t)\ dr} is plotted together with the simulated FPT density g(t). 

Figures 4-^7 show the plots of the simulated FPT density g(t) for the periodic boundary S(t) = 
So + B sin(2 7rt/Q) and of the function R[Z(t)] cxp{- /„* R[Z(t)} dr}, with R[Z(t)] given in (5.4) for 
various choices of parameters So, B and Q. Figure 4 refers to the case So = 2, B = 0.1 and Q = 3, 
Figure 5 to So = 2, B = 0.5 and Q = 3, Figure 6 to 5 = 2, S = 1 and Q = 3 and Figure 7 to So = 2.5, 
B = 0.1 and Q = 3. Note that already from rather small times, g(t) is susceptible of an excellent 
non-homogeneous exponential approximation. 

6 Concluding Remarks 

Gaussian processes play an important role in numerous fields. Although many of their properties have 
been deeply analyzed, very little exists in the literature concerning the first-passage-time probability 
density function in the presence of constant or time- varying boundaries, which would be suitable to make 
predictions on a variety of systems evolving in the presence of some critical regions of their state-space. 

In the present paper, starting from some existing contributions to this problem area, for a class of 
stationary Gaussian processes it has been proved that a non-homogeneous exponential approximation 
holds for the first passage time probability density function in the presence of boundaries that either 
possess a horizontal asymptote or are asymptotically periodic. Furthermore, for a stationary Gaussian 
process with zero mean and with damped oscillatory covariance originating at xo = 0, extensive simula- 
tions have indicated that the FPT pdf g{t) is susceptible of an excellent non-homogeneous exponential 
approximation for cither constant or periodic boundaries, even though these are not very distant from 
the initial value of the process. 

We trust that such results may prove useful for the description of the time evolution of systems 
characterized, for instance, by relaxation times much smaller than the mean observation times, thus 
making particularly appropriate and effective the asymptotic approximation obtained in the foregoing. 
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Figure 5: Same as in Figure 4 for the boundary S(t) =2 + 0.5 sin(2 7ri/3). 



A Appendix I 

Here we prove relations (3.15), (3.16), (3.17) and (3.18) of Theorem 3.1. 
Recalling the first of (2.7) and making use of (3.1) we have: 

lim ^(00= lim [e(t?i)-so7(0i)]=O (t = 1,2, . . . ,n), 

Do- >+oo 6o^+OC 

where the zero limit follows from (3.2), (3.7), (3.8) and (3.13). This proves (3.15). 
By virtue of (3.7) and (3.8), and recalling (3.13), from (2.5) one has 



1, i = j = 0, ...,n 

lim Xi+ij+i (t?i,...,i? n ) = { -7(0), i = j = n+l,...,2n 

0, i + j- 



(A.l) 



Hence, as So — > +oo the matrix A 2 „+i . . . , #„) becomes diagonal with the first n + 1 elements equal 
to unity and the last n elements equal to —7(0). Therefore, (3.16) holds. 
To prove (3.17), we first notice that 



lim (t?i, . . . ,i?„) 



[-7(0)]", i=j = 0,...,n 
[-7(0)]" i=j = n + l,...,2n 
0, * ^ J- 



(A.2) 



Therefore, by making use of (3.16) and (A.2), from the second of (2.7) one obtains 



lim ff(i?i,...,i? n ;£i,...,£„) 



lim exp, 

\ 2|A 2n+ i(i?i,...,i?„)| 



n s 
^ 'i+ntlj+ii+l( , 'lr-ii)^i 



exp 



{ 2[- 7 (0)]§^ 2 }' 
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Figure 6: Same as in Figure 4 for the boundary S(t) = 2 + sin(2 7ri/3). 



so that (3.17) follows. 

Finally, we prove (3.18). Recalling the last of (2.7) and (3.4), one obtains 

1 (27r) n 

[R(S )] n K (^---^n^l,...,U= { _ m]n/ 2 



x exp-( 



a2 



2 2|A 2n+ i(0i,...,0 n ) 



> • • • ,V n ) — 

i,3 = l 



So 



+ l i+l,j+n+l k+n+l,j+l(&l, ■ ■ ■ ,&n) g2 & 



i.j = l 



We note that from the first of (2.7) and from (3.1), for i = 1, 2, . . . , n one has 



urn = hm mzpm = lim 



So + Q(0i) 7(l?i) 

- x o — 



S 



(A.3) 



(A.4) 



where the unit value follows from (3.2) , (3.7), (3.8), (3.9) and (3.13). Taking the limit as So — > +oo in 
(A.3), and making use of (3.16), (A. 2) and (A.4), we finally obtain (3.18). 



B Appendix II 

Here we prove relations (4.31), (4.32) and (4.33) of Theorem 4.1. 

By virtue of (3.8), (4.27) and (iv) of Proposition 4.3, relation (A.l) again follows from (2.5). Hence, 
as So — > +oo the matrix A2„+i (f?i, . . . , i? n ) becomes diagonal with the first n+ 1 elements equal to unity 
and the last n elements equal to —7(0), so that one immediately obtains (4.31). 

The proof of (4.32) is analogous to the proof of (3.17) of Theorem 3.1, taking in account (4.27) and 
(iv) of Proposition 4.3. 

Finally, we prove (4.33). Recalling the last of (2.7), one has 

expj^^o + Z^)] 2 ! «■(<?!,... ,tf n ^i,...^ n ) (B.l) 
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Figure 7: Same as in Figure 4 for the boundary S(t) = 2.5 + 0.1 sin(2 7ri/3). 



Moreover, 



= exp b§ 



S 



o2 
°0 



2 |A 2n+ i (tfi,..., 0„)| 



' ' ' *o + Z(^) S S + Z&) 



So 



^) S + 



S + Z S 2 



,. S + Z(0i) 
hm — 1 — - = 1 

So— >+cx> Oq 



(i = 1,2,.. . ,n), 



(B.2) 



since is a bounded function independent of So. Furthermore, we note that from the first of (2.7) 
and from (4.1), for i = 1, 2, . . . , n one obtains 



lim 



lim 



So^+oc 5 + Z(di) So^+oc S + Z(l?j) 

5 + 



lim 



X 



So + Zfa) So + Zfa) 



1, 



(B.3) 



where the last equality follows by using condition (raj of Proposition 4.3 and by recalling (3.8), (4.2) and 
(4.22). Taking the limit as S -> +oo in (B.l), and making use of (3.16), (A.2), (B.2) and (B.3), one 
thus obtains (4.33). 
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